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ABSTRACT 



We consider the model of quantum dots interacting with coherent radi- 
ation when the relaxation processes may be neglected. The system under 
investigation consists of two discrete energy levels of the quantum dots in 
the presence of strong electron-electron Coulomb interaction and the transi- 
tions between these levels in response to electromagnetic radiation. By using 
the suitable generalisation of the Hubbard model the system of equations 
describing the evolution of the state of quantum dot was derived. 
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1 Introduction 



During the past years, a considerable attention has been paid to the growth 
and characterisation of quantum dots |], |, |], S H; Hi- The strong interest 
in these systems is readily understood in view of their obvious potential for 
future high-speed electronic devices as well as for their intriguing physical 
properties. The (semiconductor) quantum dot (QD) is system where elec- 
trons are confined in all the three dimensions. The famous Hubbard model 
provides a description for such a system || |j, 

The majority of the investigations are devoted to the transport properties 
of the coupled quantum dot system. At the same time the optical properties 
of low-dimension objects also attract the attention. For example, the possibil- 
ity of the existence of the trapped states and coherent population transfer in 
semiconductor quantum walls has been demonstrated theoretically 0, |, [10 



A solid-state implementation of stimulated Raman adiabatic passage in two 



coupled semiconductor quantum dots was proposed in [|Tl|]. The radiative 
effects of a QD array in the presence of a static magnetic field were investi- 
gated in [|12|]; with the peculiarities of the second harmonic generation in a 
two-dimensional array of QDs |13| dealt with. 

In the theory of coherent optical transient phenomena the two-level atom 
model plays an important role |fL4| . A number of phenomena (as photon 
echo, optical nutation, free- induction decay and so on) can be described in 
the framework of this model. It is attracted to develop the simplest model 
to consider the interaction of coherent radiation with QD, which is a gen- 



eralisation of the two-level atom model. In [13] the QD with only two one- 
particle energy states was proposed. There the interaction of electrons with 
scalar field of the ultra-short electromagnetic pulse (USP) was considered. If 
one takes both polarisation of the radiation and spin states of electrons into 
account, then the more complete model results. The goal of this paper is to 
derive the complete system of equations describing the two-level QD states 
evaluating under USP action. The term USP means that all relaxation pro- 
cesses in an ensemble of QD are neglected. 

The rest of the paper is organised as follows. In Sec. 2 the system of 
the equation describing the interaction of the QD with the USP is derived, 
Sec. 3 demonstrates an example of the steady-state solution of this system 
and summary is given in Conclusion. 
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2 The constituent model 



Let us consider the model of quantum dot, which has only two one-particle 
states | a) and \b) with corresponding energy levels e a and s^. We suppose 
that state \a) lies in low energy band and state \b) belongs to the higher 
energy band of the bulk material. 

The total Hubbard-type Hamiltonian of the model is written as 

H = e a %> + U aJ2 + Ta Yl PJ+i A + ^Vim) + 

in j m 

+ (1) 

+ ^ (Via+?,u + V 2 a+ 6 jU + o iA + V?6+ a iT ) + U ab ^ n^iV^. 



In the above Hamiltonian operator at (fcL) create an electron with spin /i 
in the discrete level of the one-particle states \a) (\b)) of the jth dot. The 
parameters U a , U b and U ab correspond to the electronic repulsion on the 
two states. T a and TJ, are the tunnelling matrix element between the single 
particle states. The operators n JjU = al/J^ and Nj^ = b^ ' bj^ are the particle 
number ones in the upper and lower states. Interaction of the electrons with 
electromagnetic field describes by the matrix elements V\^- The summation 
index j runs over all QD and Greek index \i runs over \ and f. 

The operators at and obey the following anti-commutation relations 

Qio-Qj/j, ~)~ [flier ^afi^iji bi a bj^ -\- b^^){ G ^aii^iji 

a ia a jfl + a jfJ a ia = 0, b^b^ + b^b^ = 0, (2) 

Oiabjn + b jtl a ia = 0, a i( jb^ + b^ ia = 0. 

In the Heisenberg pictures, operators a JM and & J/Lt obey the usual Heisenberg 
equation of motion. By using the Hamiltonian ([I]) and the relations (0) we 
can write these equations as: 

8^ 
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+ U ab a ja ^ Nja + Vi5 a J) n + V 2 5^bji, (3) 

p- 

+ U ab b ja ^ ^3° + Vi^l^tf + V* Safin- ( 4 ) 

In order to consider the one-particle energy levels population dynamics it is 
instructive to write the equations for operators n J(T and Nj a : 

d 

ih ffiNja = T b (b^bj+la -bj +1 abja^ + T b \b}j>j-l ff ~ &j"_ lo V) ~ 

- 6 al (VlS}^ - Vfi+aji) - 5 ctT (V 2 a+6 jU - V?6+ a JT ) . (6) 

One can see from these equations that population of the states \a) and \b) 
varies with time due to both tunnel processes and electromagnetic interac- 
tion determined by the transition operators Py = btfiji and P 2 j = &^%'f 
The equations (|3|) and (H) allow to write the system of equations for these 
operators: 

d - 

- (U b - E/o6)% + r a 6+ T + a i+u ) - T 6 (6+_ 1T + bf+^a^, (7) 

- {U a ~ U ab )N n P 2j + T a b+ (Vit + Vit) - T M-n + %-i|)%T- ( 8 ) 

The values of the parameters e , e&, £/ a , U b , and C/ a 6, for typical GaAs struc- 
tures used in experiments are e a ~ e b ~ 0.05 meV, C/ a ~ U b ~ 1 meV, and 
the tunnelling energy T a and T& are roughly of the same order as the level 
spacing s a . The behaviour of such a system will be crucially influenced by 
electron correlation effects because here we have a situation where U a /T a ^> 1 
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(a = a, b ) . It should be pointed out that when the dots are not separated 
far enough, the strongly localised electrons on one dot produce a significant 
potential on adjacent dots. In this paper we suppose that the distances be- 
tween QDs are so long that the tunnelling between neighbouring quantum 
dots can be neglected. It is one site approximation. The ensemble of such 
dots is similar to gas of the resonant atoms or impurities in glass. It is useful 
model being considered in the resonant optics. Henceforth the site index 
j of the operators can be omitted. Thus, we get the simpler model of the 
low-density ensemble of QDs, which is described by the following system of 
the equations 

iH^-Pi = -AeP 1 + AUR 2 P X - V X R U (9) 
at 

iH^-P 2 = -AeP 2 + AUR X P 2 - V 2 R 2 , (10) 
at 

ih-R, = 2 (^P+ - V*P X ) , (11) 

ih^R 2 = 2 (y 2 P+ - v;P 2 ) , (12) 

where Ae = (eb — e a ), AU = (U a — U a b) = (Ub — U a b). Hereafter, we suppose 
that U a = Ub 7^ U a b and define the operators Ri — — R 2 = — JVj_. 
Furthermore, in the equations @ and (]ToD the productions of the operators 
Si = R 2 P\ and S 2 = R\P 2 are appeared. By using the equations (§ - [L2D we 
can obtain 

d ~ f d ^ \ ^ ~ / d 



-AeSi + AUR 2 R 2 P 1 - V 1 R 2 R 1 + 2 (v 2 P+ - V 2 *P 2 ^j P u (13) 



ih—s 2 = ih—R 1 )p 2 + r 1 [ ih—p 2 



dt ' V dt 7 * ' x V dt 

= -AeS 2 + AUR x k\P 2 - V 2 R X R 2 + 2 (v x P+ - V x *P^j P 2 . (14) 

Here we get new operator products. Rather than continue this process, the 
algebraic properties of the operators P x and P 2 will be considered. The 
following relations result from the definition of P x and P 2 . 

P\P 2 = P 2 Pi, P x P 2 = P 2 P X , RiR 2 = R2R1, P\Ri = R2P1, R\P% = PiR\i 



5 



R 1 P 1 = -P u N ] P 1 = P u P 1 N r = 0, P 1 R 1 = P u P 1 n l = P u = 0, 

R,P+ = P+, P+N, = P+, P+nt = 0, P+R, = -P+, = P+, 

R 2 P 2 = -P 2 , = P 2 , P 2 N l = 0, P 2 P 2 = P 2 , P 2 n T = P 2 , n T P 2 = 0, 

P 2 P+ = P+, P+N x = P+, P+n T = 0, P+P 2 = -P+, ra T P 2 + = P+. 
These equalities lead to commutation relations 

P^-RiPi = 2Pi, Pi^i + ^iPi = 0, 

P 2 R 2 — R 2 P 2 = 2P 2 , P2R2 + R2P2 — 0, 

P+P1-P1P+ = -2P+, P+Px + PxP+ = 0, (15) 

P 2 + P 2 — P 2 P 2 + = — 2P 2 + , P 2 R 2 + R2P2 = 

There are more complete expressions 

P+Pi = th(1-% P 1 P+ = iV T (l-n i ), (16) 
P+P 2 = n T (l-% P 2 P 2 + = iV i (l-n T ), (17) 

which cause the following commutation relations 

P+Pi - PLP+ = (^ - iV T ) = Pi, P+Px + PLP+ = (18) 

P 2 + P 2 -P 2 P+ = (n T -iV i ) = P 2 , P 2 + P 2 + P 2 P 2 + = Rl (19) 

where the properties n^. = n CT , TVj = TV^ have been used. 

Let us denote the energy states of the QD by vectors \n a ^n a i; n^n^ >, 
where n aa (n^) is number of electrons with spin in state \a) (\b)). It should 
be note that the state of the quantum dot in one site approximation evaluates 
in the Hilbert space (& generated by four basis vectors 

1 1, 1; 0, 0, ) = |1> , |1, 0; 0, 1, ) = |2) , |0, 1; 1,0,) = |3> , |0, 0; 1, 1, ) = |4> . 

(20) 

By direct calculations in space (55 we can obtain the following representations 
for the operators R\ and P 2 . They are Pi = diag(+l, +1, — 1, — 1) and 
Pi = diag(+l, —1, +1, —1). Thus we have R\ = R\ = 1. Furthermore, one 
can obtain the following expressions 

n l = P+P 1 , n T = P+P 2 , jV T = P 1 P+ N l = P 2 P 2 + - (21) 
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By taking into account these results the reduced commutation relations can 
be written as 

P l + P l -P l P l + = R h P+P, + P Z P+=1, 
P^-RtPi = 2P h PtR + RP^O, 
P+R-RP+ = -2P+, P+P, + P z P+ = 0, (22) 

Pi Pm ~ Pm,Pl + = 0, PlP m ~ P m P\ = 0, 
P\Rm ~ Rm.Pl — 0, R\R m ~ Rm,Rl = 0. 

where l,m — 1,2. Thus the constriction of the considered algebra of the 
operators jp,P m j on the space (25 generates the Lie-algebra su(2) x su(2). 

Taking the commutation relations Q2"2] ) into account we can find R\ = 
R\ = 1. Let us introduce new operators 

N 3 = R t R 2 = R 2 R t , W = P+P 2 = P 2 P+, K = P 1 P 2 = P 2 P 1 . (23) 

These operators are governed by the equations 

ih-N 3 = 2 (y 2 P+ - V*P 2 ) Rt + 2R 2 (ViPt - K*Pi) (24) 

ih—W = (-A £ P+ + AUR 2 P+ + V*R^ P 2 + 
+P+ (-AeP 2 + AUR X P 2 + V 2 R 2 ^j = 
= AU (p+PiP 2 - P 2 Pi + P 2 ) + ^*PiP 2 - V 2 P+R 2 , (25) 

ih—K = -2AeP 1 P 2 + AU (p 2 PiP 2 + PiPiP 2 ) - ViPiP 2 - V 2 PiP 2 . (26) 

By using the commutation relations one can find that 

P 1 + PiP 2 = — P 1 + P 2 , R 2 P^P 2 = P^R 2 P 2 = — P x + P 2 , 

R 2 P 1 P 2 = P l R 2 P 2 = -P 2 Pi, P1R1P2 = P1P2 = P2P1. 

Hence, the first term in equation ( p5|) and the second term in equation ( |26| ) 
vanishe. By combining the resulting equations we can represent completed 
system of the equations for relative operators: 

d ~ ~ ~ 

ih— Pi = -AePi + AC/Si - KiPi, 
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ih^-Sx = -AeSx + AUP 1 - ViN 3 + 2 V 2 W + - 2V 2 K, (27) 
at 

h^-S 2 = -AeS 2 + AUP 2 - V 2 N 3 + 2V X W - 2VSK, 
at 

ih^-W = V?R X P 2 - V 2 P+R 2 , 

ih^-K = -2 AeK - (ViS 2 + V 2 Si 
at V 

Instead of solving for the full system of operators' equations fl2~T|) one 
may find the classical (c-number's) equations describing the evolution of the 
QD variables. These equations play a similar role in the theory of coherent 
responses of the optically excited single QD much as the Bloch equations 
are used in the case of two-level atoms. To obtain the desired equations 
the operators in equations (p7|) should be substituted for them expectation 

values, i.e. A — > f° r an Y operator A. It is convenient to introduce the 

following variables and parameters: 

2(^1,2) =Pi,2, 2(s 1 ^ = s 1 , 2i 2^ = u, 2(w} = w, 



R l2 J = n 1;2 , <^N 3 J = n 3 , 

V\,2 — ^_Rei,2, = —As /2hjjR, A = —AU/2Hujr, t = 2u>Rt, 

where lor = cLi 2 Aq is picked Raby frequency, d\ 2 is a matrix element of the 
dipole-moment operator, ei^ = Ei t2 /A are the normalised electric fields 
of the electromagnetic waves with different polarisations. In terms of these 
variables we have following system of equations: 

1^— = iZpi - Asi - eini, 
or 

1— — = Up 2 - As 2 - e 2 n 2 , 
or 
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dsi 

i—- = Qsi — Api — ein 3 + 2e 2 w* — 2enU, 
or 

= Qs 2 - Ap 2 - e 2 n 3 + 2e x w - 2ejw, 



Or 



i-q^ = g ( ei ^i _ e i^) ' ( 28 ) 
■ 9n 2 1 / * * x 

% -q^ = ^ ( e 2?2 - e 2P2) , 

,dn 3 1 1 

^-^T = 2 ( eiS i ~ e i Sl ) + 2 ^ 2S2 ~ ' 

*-jjf = g ( e i s 2-e 2 s 1 ), 

9w 1 . N 

= 2i2w - - (eisi + e 2 s 2 ) . 
<7r 2 

It is worthy of note that in these equations we take not the slowly varying 
envelope and phase approximation (SVEPA), hence ei j2 = e* 2 - However, 
SVEPA is of frequent occurrence in resonant non-linear optics. To employ the 
equations (p5[) in framework of the SVEPA the variables ei j2 should be read 
as the slowly varying complex envelopes of the ultra-short electromagnetic 
pulses and parameter Q = —Ae/2hjjjR in fl28| ) should be substituted for 
Q = —(As — TujJq) /2hujR. Here ujq is the carry wave frequency. 



3 Illustrative example 

The ultra-short electromagnetic pulse propagation can be considered in frame- 
work of the equations fl2"8|) and reduced Maxwell equations 

^ = i <Pi> , fr = i <P 2 > , (29) 

where ( = z/L ab is normalised space co-ordinate, L ab is absorption length. 
As a simple example of solving the system of equations (p8[) and (|29|), we 
consider the solution that describes the propagation of a steady state solitary 
wave, i.e., the USP. Let all variables depend only on rj = r — (/a, where a is a 
parameter. It leads to system of ordinary differential equations. The solution 
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of the resulting system can be obtained if to consider the electromagnetic 
waves to be a circular polarised ones. In this case we found 

ei, 2 0,C) = (2/r 1)2 )sech[(r-C/ai,2)/'7'i,2]exp(-zOT), e 2 ,i(r,C) =0, 

(30) 

where T\p is re-normalised USP duration. These solutions represent the 
simplest 27r-pulse by McCall and Hahn [T(J propagating in the low-density 
ensemble of quantum dots. The numerical simulation of the collision of two 
steady states with different circular polarisation shows an inelastic interaction 
of these USPs, i.e., the collision between pulses results in the generation of the 
radiation from these pulses. That reflects a strong electron-electron Coulomb 
interaction. If parameter A in equations ( p8|) set zero, then the different 
polarised steady state pulses ( PD| ) are interacted as solitons (or 27r-pulse). 



4 Conclusion 

We have introduced the simplest model of the quantum dots interacting with 
electromagnetic wave. That leads to new system of equations, which can be 
considered as generalisation of the famous Bloch equations of the theory of 



coherent optical processes ||14j| . This system in combination with reduced 
wave equations for electromagnetic waves provide the base for description 
of the coherent optical phenomena in the low-density ensemble of quantum 
dots. The simplest families of exact analytical solutions have been found for 
moving steady state pulses. Collisions between the pulses were simulated as 
the tentative, showing that they interact inelastic, resulting in generation of 
radiation and emerging a new weak pulse. 
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